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Abstract 

In the analysis of e'/ e it has been traditional to consider the isospin-breaking effects aris- 
ing from electroweak-penguin contributions and from vr'^-r/, r/' mixing, yet additional isospin- 
violating effects exist. In particular, we study the isospin violation which arises from the 
u-d quark mass difference in the hadronization of the gluonic penguin operator, engendering 
contributions of an effective A/ = 3/2 character. Using chiral perturbation theory and the 
factorization approximation for the hadronic matrix elements, we find within a specific model 
for the low-energy constants that we can readily accommodate an increase in e'/e by a factor 
of two. 
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1 Introduction 



The recent measurement of a non-zero value of Re (e'/^) 0] establishes the existence of direct CP 
violation in K ^ mi decays and provides an important first check of the mechanism of CP violation 
in the Standard Model (SM). However, the value of Re(e'/e) which emerges from combining the 
recent KTeV and NA38 results [|I|] with the earlier NA31 and E731 results 0, yielding Re (e'/e) = 
(21.2 ± 2.8) ■ 10"^ d, exceeds the "central" SM prediction of 7.0 ■ 10"^ i, i by a factor of three. 
This compels us to scrutinize the SM predictions in further detail: here we study isospin-violating 
effects arising from the u-d quark mass difference. 

Isospin violation plays an important role in the analysis of e'/e, for the latter is predicated by the 
difference of the imaginary to real part ratios in the A J = 1/2 and A J = 3/2 K tttt amplitudes. 
The differing charges of the u and d quarks engender A/ = 3/2 electroweak penguin contributions, 
whereas tt^-t], rj' mixing, driven by the u-d quark mass difference, modifies the relative contribution of 
the A/ = 1/2 and A J = 3/2 amplitudes in a significant way. The analyses of Refs. H, 0, |^ consider 
the effect of the electromagnetic penguin operator as well as of the mixing of the neutral pion 
with both the r] [^] and the rj'. Recent analyses have incorporated electroweak penguins in detail, 
as reviewed in Ref. [|, pUf j. 

Here we focus on isospin-breaking effects in the gluonic penguin operator. This operator has 
always been described as purely A/ = 1/2 in nature, but this is only true in the limit of isospin 
symmetry. That is, although the short distance structure of the operator Qq ||11|], e.g., is manifestly 



A J = 1/2, the differing up and down quark masses effectively distinguish the interaction of gluons 
with up and down quarks, so that the {tiiiIQqIK) matrix element can possess a A/ = 3/2 component 
as well |jl2|. Alternatively, one can consider the (8^, 1r) operators of the weak chiral Lagrangian [0, 



1^ , which embraces operators such as a "hadronized" Qq. In this case one finds that quark mass 
effects in the octet operators appear at O(p^) in the weak chiral Lagrangian; this is an explicit 
realization of the isospin-violating effects we discuss. 

In the isospin-perfect limit, e'/e can be written in terms of the amplitudes Aq = A{K — > (7r7r)/=o) 
and A2 = A{K (7r7r)/=2) as 

e' UJ 



e V2|e 

where 



■e(l-f^), (1) 



ReA2 ImAo ImAs 

ui = : £ = ; \l = (2) 

ReAo ^ ReAo ujlmAo 



and u ~ 1/22 emerges from an analysis of K tttt branching ratios |jT5|, |1^]. The quantity ^ is 
driven by the gluonic penguin contribution, and a non-zero fl reflects the presence of A/ = 3/2 
contributions. 

We adopt the notation Qib to denote the contribution to fl generated by the u-d quark mass 



1 



difference[], which parallels the original discussions of tt^-t],'!]' mixing effects 0. Indeed we have 
fliB = ^r),ri' + ^p, where the quantity Qp is driven by isospin violation in the hadronization of the 
gluonic penguin and ^lr],rj' arises from tc^-i], rj' mixing. Using the isospin decomposition |^ 

A{K^^n+n-) = + (3) 
and introducing "Ap" to denote K ^ nn amplitudes induced by (8^, 1r) operators, we have 



with lmAp{K° mv) = (lmAp{K^ Tr+Tr") + lmAp{K° 7r°7r°))/2. The numerator of this 
expression vanishes in the absence of isospin violation. Note that a plurality of electromagnetic 
effects, such as final-state Coulomb rescattering in the 7i~^tt~ channel |jl8[, can also make the right- 



hand side of Eq.(Q) non-zero. We ignore isospin- violating electromagnetic effects all together, for 
they are small |18[, and focus on the phenomenological consequences of the u-d quark mass difference 



exclusively. At leading-order in chiral perturbation theory, the weak chiral Lagrangian does not 
contain quark- mass-dependent effects and only flr]+7]' is non-zero; however, as we show below. 



the weak chiral Lagrangian does possess such effects in 0{p'^). 

The paper is organized as follows. In Section 2 we study the K ^ ttit amplitudes at tree level 
in 0{p^) in the weak chiral Lagrangian. This general framework allows us to identify all the isospin 
breaking effects that can occur at next-to-leading order in chiral perturbation theory, albeit the low- 
energy constants are unknown. In Section 3 we consider the gluonic penguin operator. To estimate 
its contributions to the isospin breaking operators we identify in Section 2, we use the factorization 
approximation for the hadronic matrix elements. Within the factorization approximation, the 
terms of the 0{p^) strong chiral Lagrangian contribute to the O(p^) weak chiral Lagrangian. These 
low-energy constants are also unknown; in this case, however, we can use a resonance-saturation 



model to estimate them, to illustrate the effect. 



2 Chiral Lagrangian Analysis 

Chiral perturbation theory forms a natural framework in which to discuss isospin violation in 
the K TTTT amplitudes. The weak chiral Lagrangian is realized in terms of the unitary matrix 
U = exp(2i0//), which transforms under the chiral group SU{3)l x SU{3)r as U RUL'^ , where 
R, L are elements of 5'f/(3)/j,L respectively. The function represents the octet of pseudo-Goldstone 
^We adopt this notation for simplicity and refer the reader to Ref. for a detailed discussion. 
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bosons, where = J2a=i,...8 ^a4'a "^^^ chiral Lagrangian is constructed in terms of U and its 
derivatives, note that = iU'^ D^U , as well as in terms of the function x, which transforms as U 
under the chiral group. In the absence of external fields, x = 2BqM, where M denotes the quark 
mass matrix, M = diag(mtj, m^, m^), and the parameter Bq, proportional to the quark condensate 
(qq), has dimensions of mass. The function x encodes the isospin- violating effects of interest as 
rriu ^ nid. 

The C(p^), CP-odcl0, weak chiral Lagrangian transforming as (8^, 1r) under SU{3)l x SU{3)r 
has only one term [|l3l 

d^^ = c^{XrL'), (5) 

where C2 is a parameter of order of the Fermi constant Gp and () denotes a trace over flavor indices. 
The x-dependent term anticipated from the form of the leading-order strong chiral Lagrangian, 
proportional to {\j{x^U + f/^x)), does not appear in the computation of physical amplitudes |n . 
The term in question can be written as a total divergence in our case and thus does not 
contribute, in accord with Ref. |21[ |. 

In next-to-leading order x-dependent terms are possible. The complete weak chiral La- 

grangian has been constructed in Ref. |T^. Collecting the terms which can evince isospin violation, 
we find 

^'wliB = E^{Xrxl) + E,{Xjx+){x+)+E^{Xrx'-) 
+ E^{\jX-){X-) + E^{\7t[x+,X-]) 

+ E^,{Xr{x+,L'}) + E^,{\jL,x+L^) + E^,{XrL^){{L^,X+}) 

+ E^,{\rL'){x+)+E^,{Xrt[x-,L']) , (6) 

where x+ and X- are defined as 

x+ = x^u + u\ 

X- = i[x^\J-U\). 

We can use this Lagrangian to calculate fip from Eq. (^. Working to leading order in isospin 
breaking, so that merely terms linear in — are retained, and dropping terms suppressed by 
Ml /Ml, we find 

2^2 M2 So(md - rriu 
0.12GeV^ 



0„ — _I ill )ll 'iL[ 0E~ — IE" — AE~ — Fr — Fr — AFr — E 



Co 



2E^ - 2E3- - AE, - E^o - i?n - '^E^2 - ^fs • (7) 



^We thank G. Colangelo and J. Kambor for their generous assistence in rectifying the notational errors of our 
original manuscript. 

■^We choose to drop O13 as the one of the O10-14 which is not independent because it does not contribute to 
Eq. (H). Also, O32-34 are not included because they are related to Oi_5 by equations of motion. 
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The terms proportional to E^, E2 , and can potentially generate tadpole contributions, as the 
eigenstates of the weak interaction are not those of mass. Consequently, we take care to remove 
possible tadpole contributions via the construction of Ref. |T^. Note that to 0{md — rriu) and 
to the order in the momentum expansion to which we work, it suffices to use Eq. (^) to compute 
Im Ap{K^ — » vrvr). Thus merely C2 appears in the denominator of Eq. (|^. The numerical value 
given reflects the use of Ref. ||2^ and of the relation 



Boirud - m„) = M^o - + - M^o (8) 



which follows from the leading-order strong chiral Lagrangian, in concert with Dashen's theorem [^, 
2^. The utility of Eq. (|^) is limited, for the E~ coefficients are unknown. However, power counting 
in chiral perturbation theory suggests that each of the constants E^ is suppressed by C(A^gg) with 
respect to . Thus the numerical prefactor in the last line of Eq. ought determine the "natural" 
size of Qp — it is of order 0.1. Remarkably, these effects are of comparable numerical size to the 
value of flf] in O(p^) P, 01) so that the terms found in Eq. (J^) merit further study. 

For reference, it is useful to summarize the results of Refs. ^ for flr]+r]' in C^(p^) and then to 
proceed to enumerate all possible isospin-violating contributions in 0{p^) , irrespective of whether 
we term them or "i7p". In leading-order chiral perturbation theory, ir^ — t] mixing is the 

only rrid 7^ rriu effect to impact the K —>■ mr amplitudes. Diagonalizing the neutral, non-strange 
meson states of the strong chiral Lagrangian in 0{p'^) , noting 

4'^ = f f(^"^.) + (x+)) , (9) 



(10) 



4 

yields the physical vr^ state in terms of the octet fields vr'' and t] 



TT = vr^ + — — T] + 0{md - rriu] 

\ J phys 4 \ rris — m / 



where m = (m^^ + mu)/2. Using Eq. (|^) to compute lmAp{K 'n-^ri)/lmAp{K 7r°7r°) yields 
finally ^ 

^ 1 / 171.1 — m,, \ , , 

noting {ms—rh) / {md—mu) = 40.8±3.2 The approximate equality of this result to the numerical 
coefficient of Eq. (^ follows as 4:{M'^ — M^) ^ O(lGeV^). Consequently, it is also important to 
evaluate the impact of tt^-i] mixing in 0{p'^) on the K ^ Tin amplitudes. This has already been 
done to some extent, for the usual analysis 0, 0] deviates from strict chiral perturbation theory in 
that an explicit r]' degree of freedom appears as well, leading to both 7r° — rj' mixing and t] — r]' 
mixing. For comparison, explicit study of ir^-r] mixing in 0{p'^) , noting [24 
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shows that it is sensitive to the low-energy constant L7 EM. In an effective Lagrangian which 



includes the rj' degree of freedom via the nonet symmetry of a large Nc approach, taking the limit 
of small momenta and <C M^, yields an interaction of the form associated with L7 



Moreover, the rj' contribution numerically saturates the value of L7 found phenomenologically [0, 
The presence of the t]' thus apes higher-order effects in the strong chiral Lagrangian. Including 
the rj' as per the usual analysis [§, yield^ 

= fiJ (cos^ - V2sine)2 + ^ ^(sin^ + y2cose' 

^ 2.4 fi^ ^ 0.31 (13) 

where we use 9 = —22° for the rj — rj' mixing angle 0, ^ . The effect of the rj' is no smaller than 
that of the rj; this is consistent with the comparison of Eq. (^ with Eq. (|TT|).Q There are thus a 
plurality of effects which are important in O(p^) . Let us enumerate the possible isosp in- violating 
effects which occur in 0{md — m^) and 0{p'^) : 

i) Isosp in breaking in the 0{p^) mass term of Eq. (^, including tt^-t] mixing, acting in concert 
with the O(p^) weak chiral Lagrangian, Eq. (|^), computed to one-loop order. 

ii) 71^-7] mixing, realized from the 0{p'^) mass term of Eq. (P), combined with the isospin- 
conserving vertices of the 0{p'^) weak chiral Lagrangian. 



iii) Next-to-leading order n^-rj mixing as per the strong chiral Lagrangian in 0(p^), Eq. (|T2|), 
combined with the leading-order weak vertex from Eq. (^. The vr"-//' mixing effects of the 
usual analysis are an example of this type. 

iv) Isospin violation in the vertices of the 0{p'^) weak chiral Lagrangian, Eq. (^. This is realized 
as Eq. (1^ and serves as our focus here, for it contains the qualitatively new effects we argue. 

We wish to focus on the contribution of iv), yet we cannot avoid considering that of ii), for the 
low-energy constants Ei of Eq. (|^) potentially enter here as well. Considering exclusively the terms 
of Eq.(^ we find the contribution of ii) to be: 

12GeV^ / 1 3 \ 

-2(^3- + - E^) + -(Efo - ^n) - '^E^2 + E^A + TT^fs , (14) 



"^Note that using the ir^-rj, rj' mixing formulas resulting from the exact diagonalization of a chiral Lagrangian 
based on nonet symmetry in OijP') and large Nc p8| yields ^n+ri' — 1-7 £7^ ~ 0.22. For comparison, note ri,,+,,' = 
0.25 ± 0.02 from the recent analysis of Ref. 

^Large Nc arguments suggest that L-j could dominate the low-energy constants in 0(jy^) ||2^, yet this is phe- 
nomenologically not the case 0, |2^ . 
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so that there is no manifest cancellation with the terms of Eq. (|^). Moreover, we find in Section 3 
that the contributions of Eq. (11) are numerically smaller than those of Eq. (|^) — it is the latter 
which contains isospin-breaking effects in the hadronization of the gluonic penguin operator. To 
study these effects in detail, we must turn to the factorization approximation and estimate, as in 
the next section, the contributions of the gluonic penguin to the E~ of Eq. (|^). 



3 Factorization 

Within the context of the factorization approximation, the bosonized form of the Qq penguin 
operator appears as the product of scalar and pseudoscalar densities obtained from the strong 
chiral Lagrangian. The construction relevant to tttt decay is p9l [11 



jCp = -^V*^VudC6(^-8{sLqR)iqRdL)^ + h.c. 

• ^^v:,VudCe32B',^^ + h.c., (15) 



Ci2 



where q{L,B) = (1 =F 75)q'/2 and Cq is defined as in Ref. ||3^. Using the 0{p'^) strong chiral La- 
grangian p4|, Eq. (|12D, one finds 

C2 = ^CKJmCe (iQB'ofnL,) • (16) 
Equation ([T5|) also yields a term of the form {\t{x+))j proportional to Lg, yet this is merely the 



weak mass term discussed earlier — it does not contribute here ||20| , |2T[| . It is well-known that the 
contribution of the CP-even analogue of Eq. ([I6|) does not suffice to reproduce the phenomenological 
value of C2, its associated low-energy constant |2^, where we note ImCg — > ReCg in Eq. (|16D yields 
C2 ~^ C2. Equation (p!6D is useful nevertheless, for it serves to normalize the isosp in- violating 
constants induced by the Qq operator. 

The O(p^) strong Lagrangian, Eq. (|T^), as per Eq. (0), also yields contributions to certain of 
the coefficients enumerated in Eq. (^, as well as to other operators of the 0{p'^) weak chiral 
Lagrangian. The non-zero contributions to are 



^10 - ^15 - f2 ' ^13-— FTF ' ^14 - f2 ■ y^') 

In In -^5 In 

Unfortunately, however, this approach does not yield a complete estimate of the coefficients of 
the O(p^) weak chiral Lagrangian in the factorization approximation. The bosonization of Qq, as 
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defined in Eq. (|1^), also contributes to the 0{p'^) weak chiral Lagrangian through the 0{p^) strong 
chiral Lagrangian. Although the latter has been constructed [31, 19], its coefficients are not known, 
and we must turn to a model to proceed. 

The use of resonance saturation allows us to estimate some of the coefficients in the 0{p^) strong 
chiral Lagrangian. This has been done for the Lj constants that appear in the 0{p^) strong La- 
grangian [^|. In particular, the form of the terms needed in Eq. suggest that scalar and 
pseudo-scalar resonances might be dominant. This is true for the coefficients of the 0{p'^) strong 
Lagrangian that appear in Eq. ([iTl). The constant Lj is saturated by the rj' [^, and it is reason- 
able to assume, in an analogous manner, that the constants L5 and Lg are saturated by the scalar 
resonances . Indeed, Ref . ||26| inverts this argument and uses the phenomenological values of L5 



and Lg to fix the couplings of the scalar resonances to the pseudoscalar octet of vr's and 77's. 
As a model for the needed 0{p^) counterterms, we propose the Lagrangian 

Cs = \{D^SD,S - Mis') + c,{S}SiL') + c„(e^5^x+) + ^{i^S'ix^) (18) 

for the scalar meson octet, where t/ = The ffist three terms of this Lagrangian are explicitly 
considered in Ref. In the limit of momenta such that -C M|, the scalar octet no longer 



plays a dynamical role and is thus "integrated out," yielding |^ 



2 

- Mi" ' - 2Ml- 

Using a scalar mass of Ms = 0.983 GeV and assuming that this contribution saturates Ll^^Mp), 
one finds = 0.042 GeV, q = 0.032 GeV |2l. 



The last term in Eq. ([Iq ) has been recently considered in Ref. |T9[. This term breaks the mass 
degeneracy of the states in the scalar octet, splitting the Kq{14:30) from the ao(980), for example. 
This term also generates some of the 0{p^) strong operators of interest to us. Integrating out the 
scalar octet one finds two terms proportional to dm in 0{p^) 



~ 2M| + ~Mr ' ^ ^ 



which contribute to the scalar densities in the bosonization of Qq. The new contributions are 

= ^ _4.8 GeV-2 ; ^ = ^ -2.4 GeV'^ . (21) 



C2 2MIL5 C2 MIL, 

Within the context of our model, other contributions of the 0{p^) strong Lagrangian to the constants 
Ei are assumed to be identically zero. For the numerical estimates we ffi dm to the K^-ao mass 
difference, 

Ml* - Mj, 
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and use L^fi = Mp) of Ref. Note that the terms of Eq. (|21| ) are numerically much larger than 
those of Eq. ([TtI) — indeed, they dominate VL-p. If we use the values of L5 7 g(/i = Mp) from Ref. p6[| : 
L5 = 1.4 X 10-^ L7 = -0.4 X 10-3, and Lg = 0.9 x lO^^, Eq. (O) yields 



^ = 0.13 GeV-2 ; ^ = -0.24 GeV' ; ^ = 0.21 GeV^' , (23) 
C2 C2 C2 

where = 93 MeV. In view of the dominance of the terms computed from the 0{p^) coefficients, it 
is important to compare the relative size of the 0{p^) and 0{p^) coefficients induced by the scalar 
resonance. To illustrate, let us consider the ratio of the coefficient of the first term of Eq. (pOf) , 
calling it to the coefficient Lg: 

This ratio is large, but not inconsistent with dimensional analysis. Nevertheless, it may be naive 
to associate the Kq-qq mass difference with flavor-symmetry breaking as in Eq. (pISl). For example, 
quark model studies suggest that the ao(980) may well be a KK molecule [^, If we use the 
predicted lowest-lying isovector and strange scalar states of Ref. [Q, yielding masses of 1.09 and 
1.24 GeV, respectively, we flnd, rather, that dm ~ —0.76 and that /^i/Lg ~ —0.79 GeV"^ — this is 
also consistent with dimensional analysis. 

We can now proceed to estimate the value of f2p using our estimated low-energy constants. Were 
we merely to use the numbers of Eq. (^Tj) and = —2.4 we would obtain 

Qp = (0.12 GeV^) ^"^^"^(3^"^ - ""d) _ _Q g5 (25) 

s 5 

This unexpectedly large result is driven by the value of dm found in Eq. (^): using dm ~ —0.76 
yields Qp ~ —0.28. The sign of dm and thus of Qp in our picture is the consequence of the mass of 
lowest-lying strange scalar being greater than that of the lowest-lying isovector scalar. Collecting 
the contributions of Eq. (|2l|) and Eq. (^3|) yields 

= -0.79 (-0.21) (26) 

for dm = -2.4 (-0.76). Note that Eq. (|T|) and L^(/i = Mp) = 0.2 • 10"=^ [||] yields n^l^, = 
-0.12 (-0.03), so that writing fi^ = + fil,'?^ yields = -0.91 (-0.24). For reference, the 
value of Eq. (^ used in the "central value" of e'/e in Ref. is Q^ri+^n' = 0.25±0.05, whereas that used 
in Ref. is ^lr]+v' ~ 0-25 ± 0.10. The changes in fli^ found in O(p^) impact e'/e in a signiflcant 
manner. Using the simple formula of Eq. (1.7) in Ref. ^ shows that under Qib = 0.25 —>■ —0.25 
the value of e'/e increases by a factor of 2.2 . Thus a very small or negative value of Qib generates 
an increase in e'/e with respect to the usual value cited [Q. It is particularly noteworthy that the 
range in our estimates of fijg exceed the central value of Qn+ri' used in Refs. |Q, |T^. The detailed 
numerical results we find do rely on a simple model; nevertheless, a substantial increase in the error 
associated with Qib, Eq. (^, is in order. 
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4 Conclusions 



We have shown that there are isospin-breaking effects in e' which have not been previously con- 
sidered. Specifically, we have examined the role of isospin violation in the matrix elements of the 
gluonic penguin operator within the context of chiral perturbation theory. Although the presence 
of unknown low-energy constants implies that we lack a reliable way to calculate these effects, 
we believe such limitations underscore the need for a larger uncertainty in the theoretical value 



of e'/e than currently in vogue. In particular, the recent reviews of Refs. @, use 0.25 ± 0.05 
and 0.25 it 0.10, respectively, for the value of Qib = ^r^+r^'- Our estimate of Qib from the specific 
iTT'd 7^ fTT'u effects we consider ranges from 0.1 —0.7. This range reflects a variation in e'/e of more 
than a factor of two. 
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